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Abstract. In-flight measurement of the antenna main 
beams of the Planck instruments is a crucial input to 
the data analysis pipeline. We study the main beam re- 
construction achievable through external planets using a 
flight simulator to model their observation. We restrict 
our analysis to the 30 GHz LFI channel but the method 
can be easily extended to higher frequency channels. We 
show that it is possible to fit the antenna response from 
Jupiter and Saturn to obtain an accurate, robust, simple 
and fast reconstruction of the main beam properties under 
very general conditions, independently of the calibration 
accuracy. In addition, we find that a bivariate Gaussian 
approximation of the main beam shapes represents a sig- 
nificant improvement with respect to a symmetric repre- 
sentation. We also show that it is possible to combine the 
detection of the planet's transit and Planck's very accu- 
rate in-flight calibration to measure the planet's temper- 
ature at millimetric wavelengths with an accuracy at the 
% level. This work is based on Planck-LFI activities. 

Key words: Cosmology: cosmic microwave background - 
Space vehicles - Telescopes - Methods: data analysis - So- 
lar System: general 



1. Introduction 

The Planck Surveyor]] is the ESA space mission de- 
voted to the study of the Cosmic Microwave Background. 
Planck will have an impact on a number of scientific is- 
sues, such as the physics of the early universe, structure 
formation theory and cosmological parameters determina- 
tion (Bersanelli et al. 1996). In order to reach the neces- 
sary level of sensitivity it is important to understand sys- 
tematics and to keep them under control. In this paper we 
will focus on the behavior of the Planck Low Frequency 
Instrument (LFI, Mandolesi et al. 1998) antenna patterns. 
For simplicity, we will restrict our analysis to the 30 GHz 
LFI channel but the method we present here can be easily 
extended to higher frequency channels. 

Send offprint requests to: burigana@tesre.bo.cnr.it 
1 http:/ /astro. estec.esa.nl/SA-general/Projects/Planck/ 



The beam pattern is affected by optical distortions, 
which depend on the telescope design and on the arrange- 
ment of the various feed horns in the focal plane. These ef- 
fects degrade both angular resolution and sensitivity (e.g., 
Mandolesi et al. 2000a, b). Therefore, accurate measure- 
ment of the beam pattern is a crucial input to the data 
analysis pipeline. 

Due to their small angular size external planets pro- 
duce large signals only when seen in the main beam. As 
such, they represent a unique possibility to recover di- 
rectly from the data the in-flight behavior of the main 
beam. This point was already addressed, in the framework 
of the Planck mission, by Bersanelli et al. 1997 for the 
simple case of a Gaussian symmetric antenna response. 
We extend here this analysis to quantify our ability to 
reconstruct a more realistic, asymmetric beam pattern. 

The plan of this paper is as follows. In Sect. 2 we 
describe our main tools and assumptions. In Sect. 3 we 
discuss the quality of the main beam reconstruction and 
the implications for the treatment of Planck data. As a 
byproduct of our simulations, we also verify in Sect. |]thc 
validity of the Planck optical design to strongly suppress 
straylight contamination from internal Solar System bod- 
ies. In Sect. 5 we summarize our findings and draw our 
conclusions. 



2. Method 

In order to attack the problem of the in-flight main beam 
reconstruction, we have to: (i) describe the Planck orbit 
and scanning strategy; (ii) quantify the antenna response; 
(Hi) exploit the planet's mm emission and positions; (iv) 
simulate the Planck observations of the external plan- 
ets (basically Jupiter and Saturn) . Here we briefly discuss 
these points separately. 

2.1. Planck orbit and scanning strategy 

The selected orbit for the Planck satellite is a Lissajous 
orbit around the Lagrangian point L2 of the Sun-Earth 
system (e.g., Bersanelli et al. 1996). In the nominal oper- 
ation scheme the spacecraft spins at 1 r.p.m. around an 
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axis kept parallel to the ecliptic plane. Every hour the spin 
axis is moved by 2.5' maintaining its anti-solar direction. 
The telescope optical axis is at an angle a from the spin 
axis direction. The spin axis might precede about the anti- 
solar direction, with a period of about six months and an 
amplitude of about 10°. This spacecraft movement is of 
course over imposed to the Lissajous orbit and to the spin 
axis hourly shift. In this paper we consider values of a 
between 80° and 90°, i.e. about the value of 85° recently 
recommended by the Planck Science Team. We make 
use of the Planck flight simulator described in detail by 
Burigana et al. 1997, 1998 and Maino et al. 1999 prop- 
erly modified to model the Planck observations of the 
Solar System bodies and the spacecraft motion (see, e.g., 
Bersanelli et al. 1997). 

For what follows, it is convenient to introduce a tele- 
scope "reference frame" (hereafter rf) {xt,Ht,zt} with 
the zt axis coincident with the direction of the telescope 
line of sight (the p direction, say) and with the xt ~ Dt 
plane identifying the telescope field of view plane (we 
choose to orient xt towards the intersection of the xt — Vt 
plane with the spin axis s or, in the case a = 90°, xt \\ s). 
For the considered scanning strategies the spin axis and 
the telescope directions (s and p, respectively) are easily 
derived given the observation time, the spinning frequency 
and the boresight angle a. So it is always possible to pass 
from a chosen celestial rf to the telescope rf (and vicev- 
ersa) by a suitable Eulerian rotation of the considered rf. 

2.2. Antenna angular pattern 

The Planck High Frequency Instrument (HFI, Puget 
et al. 1998) is located at the center of the focal plane. 
The LFI feed horns surround HFI and are then substan- 
tially off-axis. For instance, with a telescope of 1.5 m class, 
the 30 GHz beams are at about ~ 5° from p. So, it is con- 
venient to define a beam rf {xb, yb, Zb} with the Zb axis 
coincident with beam axis b and with the Xb — yb plane 
slightly tilted with respect to the telescope field of view 
plane (we keep the convention of obtaining the beam rf 
from the telescope rf through a rotation of the telescope 
rf about an axis hortogonal to the plane identified by p 
and b by the angle necessary to transport p to b. 

As we know the position of each feed horn in the fo- 
cal plane, we can always pass from the telescope rf to the 
beam rf and viceversa. In fact, the flight simulator deter- 
mines for every time step the orientations in the sky of 
the telescope and beam reference frames, to compute the 
antenna response for a given line of sight. 

2.2.1. Main beam 

Recent improvements on the Planck optical design based 
on aplanatic solutions (Mandolesi et al. 2000b) show that 
the main beams are roughly elliptical, with an elliptic- 
ity ratio r ;$ 1.4 (Alcatel, private reference, PL-AS-TN- 



022). Therefore, we approximate the antenna pattern as 
an off-axis bivariate Gaussian beam. For simplicity, we will 
consider the bivariate Gaussian beam projected onto the 
field of view plane (i.e. we will consider this beam repre- 
sentation on the xt — Vt plane, and not on the Xb — yb 
beam plane). To be accurate one has to say that if the 
true beam shape is elliptical in the Xb — yb plane, it gets 
distorted by the projection on the xt — yr plane. How- 
ever, since the off-axis angle even for the 30 GHz beam is 
small, this distortion is negligible and, if anything, does 
not change the elliptical nature of the beam response. In 
addition, a realistic main beam distortion implies a devia- 
tion from the elliptical shape larger than that introduced 
by this projection. So, let (x T ,y T ) identify the projection 
of the beam centre unit vector onto the xt — Vt plane. 
Let e be the angle between the xr-axis and the princi- 
pal axis of the bivariate Gaussian. The (normalized to the 
maximum) beam response can be then expressed as: 

J = exp[-~ (TZuf^TZu] (1) 
where 

S = diag(cr^_, o~^_) contains the bivariate's beam disper- 
sions along the ellipse principal axis and 1Z is the rotation 
matrix for an angle e in the xt — Vt plane: 

I — sin e cos e J 

It is then convenient to define the beam "sigma" a = 
^/(7 + er_ and the ellipticity ratio r = <j + /<t_ . 

2.2.2. Far side lobes 

A realistic description of the antenna far side lobes has to 
rest on the accurate optical calculations of de Maagt et al. 
1998 for the Planck telescope including shields. Several 
cuts at constant azimuthal angle </> are shown in Fig. 1 as 
a function of the colatitude angle 9 from the beam axis b. 
In the "antispillover" region [about (9,<j>) ~ (90°, 180°)] 
where it is mostly important to evaluate the effects of 
straylight from the internal bodies of the Solar System 
(i.e. Sun, Earth and Moon) the pattern response drops 
down to approximately ~ lOOdB. At such a level of rejec- 
tion, according to the Planck requirements, one expects 
that Sun, Earth and Moon are completely harmless for 
the mission. We will further discuss this point in Sect. |]. 

2.3. Planet's mm emissions and positions 
2.3.1. External Planets 

Several authors reported measurements of the planets 
brightness temperature at millimeter wavelengths with 
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Fig. 1. Several cuts of the full antenna pattern, in c!B, as 
computed by de Maagt et al. 1998. The lines refer to the 
antenna response for 8 between 0° and 180° and, from the 
bottom to the top, with <f> at steps of 30° from 0° to 360° 
(each cut is vertically shifted for graphic purposes to have 
the value at the point at 9 = 0° equal to the considered 
value of (f). 



typical uncertainties of 3 -j- 5% (see, e.g., Bersanelli et al. 
1997 and references therein). The quite large uncertainties 
associated with these values prevents one from using plan- 
ets for accurate temperature calibration of the Planck 
time order data. This will be done, to better than a 1%, by 
using the diffuse signature of the CMB dipole anisotropy 
(Bersanelli et al. 1997). However, for the purpose of beam 
reconstruction it is not necessary to have a detailed knowl- 
edge of the planet emission. It only matters that the source 
is stable and sufficiently bright to be detectable even when 
the source is far from the beam axis. This requirement is 
crucial to sample the antenna beam response at different 
angles. For this reason we will consider here only Jupiter 
and Saturn, which are the brightest of the external plan- 
ets. On the basis of the published data, we will assume 
hereafter that Jupiter and Saturn have, at 30 GHz, bright- 



ness temperatures of T^ p 
respectively. 



152 K and T, 
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133 K, 



2.3.2. Inner bodies of the Solar System 

The amount of straylight contaminations from the inner 
bodies of the Solar System can be easily quantified. We 
restrict ourselves to the Sun, Earth and Moon, because of 
their high intrinsic temperatures and because of the solid 
angles subtended from L2. We will use = 6000 K, 



300 K and T^oL 



L earth — " u " iv ^ moon — 250 K as the temperatures 
associated with the millimetric thermal emission of Sun, 



Earth and Moon, respectively. 



2.4- Planck observations of the external planets 

We use the Planck flight simulator in order to model 
the transit of the planets in the Planck field of view. In 
particular, Jupiter and Saturn will be observed twice in 
about a year. The solid angle of the external planets as 
seen by Planck is very small compared to the beam size. 
Thus, the Planck observations of the Jupiter (Saturn) 
will yield 



' 30GHz 



mm 



T^ir(R/df J[j{t) - o] 



(2) 



In this equation T 30 chz is the observed Jupiter (Saturn) 
brightness temperatures @ 30 GHz; , R and d rep- 
resent the intrinsic Jupiter (Saturn) brightness tempera- 
ture, radius and distance, respectively; J is the antenna 
response and j(t) identifies the angular position of the 
planet as seen by Planck, the time dependence being 
fixed by the scanning strategy. 

3. In-flight recovery of the main beam pattern 



The Planck Time Ordered Data (TOD) are affected by 
instrumental noise. Therefore, our capability to recover 
the main beam pattern rests on the possibility to clearly 
detect a bright source (e.g., Jupiter), even when signif- 
icantly far from the beam axis. A proper description of 
the Planck-LFI instrumental noise should in principle 
include a 1// contribution (see, e.g., Bersanelli et al. 1996, 
Seiffert et al. 1997). The knee- frequency of the 1// noise 
is expected to be comparable with the spinning frequency. 
However, it has been shown that destriping algorithms can 
very efficiently remove this low frequency noise component 
even under more pessimistic conditions (see, e.g., Maino 
et al. 1999 and references therein) and return a TOD that 
we will assume, accordingly with the goals of this paper, 
white noise dominated. So, in what follows we will model 
the TOD noise component as pure white noise, with the 
Planck goal sensitivities discussed by Bersanelli et al. 
1999 (private reference, Planck Low Frequency Instru- 
ment, Instrument Science Verification Review, October 
1999, LFI Design Report). In principle, the signal fluc- 
tuations introduced by CMB and foreground anisotropies 
behave as a noise source in this context. However, since 
they can be accurately subtracted from the TOD by using 
the Planck final maps, we neglect them in what follows. 
We adopt here a simple scanning strategy with a = 90° . In 
the simulations presented in Sect. 3.1, we oversample each 
scan circle every ~ 5' (i.e. roughly 6 points per FWHM @ 
30GHz) and shift the spin axis by 5' every two hours. Af- 
ter simulating the Jupiter and Saturn transits we extract 
from the time ordered scans a few (~ 100) chunks com- 
prising the source transit. Since the source is pointlike, 
these chunks give, when displayed one after the other and 
having taken into account the small variations of planet 
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Fig. 2. Top panel: signals of different components along 
the scan circle with the maximum of Jupiter contribution 
during its first transit. Solid line: Jupiter signal; crosses: 
white noise; diamonds: 1/ / noise coupled to white noise; 
dotted line: signal from CMB and extragalactic source 
fluctuations and Galaxy emission modelled according to 
Burigana et al. 2000a and references therein. Bottom 
panel: the same as in the top panel, but considering the 
signal at the same scan position, where Jupiter signal is 
maximum, for different scan circles. According to the sim- 
ulation parameters, by multiplying the scan circle num- 
ber or the sampling number on the scan circle by 5 we 
have respectively the angular displacement between differ- 
ent scan circles and (approximately, owing to the off-axis 
beam position) the angular displacement between different 
samplings along the scan circle, expressed in arcmin. We 
consider here the case of an elliptical beam with r — 1.3, 
clearly visible in the different spread of Jupiter signal in 
the two panels. Note that the signal to noise ratio is larger 
than unit up to —(20 -j- 25) dB, i.e. at about (3 -j- 3.5)cr. 
[Signals in dB normalized to the maximum Jupiter signal 
at its first transit]. 

distance in the different samplings, a 2-D plot of the beam 
profile, ~ 8.3° x 8.3° wide. 

In Fig. 2 we show the expected signal from Jupiter as 
seen along the scan circle which crosses the source at the 
maximum (top panel) and as seen along an arc orthogonal 
to this circle (bottom panel). The signal to noise can be 
improved, wrt the case of a single receiver and transit, by 
considering that two LFI receivers are coupled to the same 
optical beam, that there are two 30 GHz beams with the 
same optical properties and that two (three) transits of 
both Jupiter and Saturn are expected for a one year (for 
a 14 -7- 15 months) mission. This obviously increases the 
signal to noise ratio by a factor of 2y2 (2v3). As a result, 
@ 30 GHz, the shape of the main beam can be recovered 
down to -(25 4- 32.5) dB, i.e. at about (3.5 4- 4)cr. 



3.1. Recovery of the main beam parameters 

We consider both a symmetric and an elliptical beam, with 
ellipticity ratio r = 1.3. The numerical values of the beam 
parameters are shown in Table 1. We use Eq.(l) to model 
the antenna response. We fit the beam shape theoretical 
parameters to the 2D plot of the beam response obtained 
as mentioned at the end of the previous section. 

The results of the fits are shown in Table 2. We fit 
also an additional parameter, r^, related to the planet 
brightness temperature and to the average distance of the 
planet from the spacecraft, (d), for the points considered 
in the fit: r k = n(R/(d)) 2 T^/ J 4?r J(t) dQ. . We recover 
the full set of parameters with very high accuracy (Buri- 
gana et al. 2000c), the x 2 /DOF being always very close 
(to better than 1%) to the unity value. It is obviously 
more efficient to recover the beam pattern parameters us- 
ing Jupiter rather than Saturn, simply because Jupiter is 
brighter. 

An interesting byproduct of the beam fitting procedure 
is the possibility of estimating the planet's antenna tem- 
perature, or more properly the product T^R 2 relevant 
for the planet's emission at the considered frequency. This 
quantity, as previously stated, is a poorly known quantity 
at LFI's frequencies. By considering Eq. ([!]) it is clear that 
the latter temperature, is related to the normaliza- 
tion of the bivariate which is in turn very well constrained 
by our fit. The ability to estimate T^ b ' then rests upon the 
overall calibration accuracy^ and on the knowledge of the 
total antenna beam integral. Calibration for Planck will 
be provided by continuous observation of the CMB dipole 
signature and its modulation introduced by the spacecraft 
motion (Bersanelli et al. 1997) and is expected to be ac- 
curate to within 1%. The total integral of the antenna 
pattern poses a more serious problem, as the contribution 
of the far side lobes goes undetected when using a celes- 
tial source. However, optical calculations (de Maagt 1998) 
show that the contribution to the antenna pattern coming 
from outside the main lobe is expected to be < (2 4 3)% 
of the total. The uncertainty on the latter figure may then 
dominate and pratically set the accuracy on the LFI es- 
timate of TW; therefore, even a poor knowledge of this 
contribution (e.g., with an accuracy of ~ 30%) allow to 
reach a ~ 1% level of accuracy in the measurement of 
T<» (Burigana et al. 2000c). 

So far no systematic effect has been included in the 
simulations. One of the main goals of the Planck mis- 
sion is to limit the contamination arising from system- 
atics at few /iK level. However, systematics affecting the 
spacecraft pointing and rotation, which are not expected 
to significantly degrade the sky temperature estimation, 
may in principle turn out to be harmful when reconstruct- 
ing the beam shape. To investigate this possibility we have 
extended our simulation to include these effects. 

2 We want to stress that the calibration of the TOD is not 
needed in order to reconstruct the other beam parameters. 
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3.1.1. Pointing uncertainty 

A potential problem for in-flight beam reconstruction is 
posed by the telescope pointing uncertainty, i.e. the lim- 
ited accuracy in determining the spacecraft's effective spin 
axis. However, it is a Planck requirement that such an 
uncertainty be less than 1'. This is why, to simulate this 
effect we introduce a pointing uncertainty drawn from uni- 
form distribution of values between and 1. Since the 
spin axis direction is not expected to significantly change 
within the same scan, the uncertainties only arising when 
the spacecraft is repointed, we keep the perturbation con- 
stant along the scan. 

The impact of pointing uncertainty in the elliptical 
beam parameter recovery for the first transit of Jupiter 
and Saturn is also shown in Table 2. Clearly, the higher 
signal to noise ratio in case of Jupiter makes the results 
more sensitive to this kind of systematics. This is evident 
from the larger degradation of the beam parameter recov- 
ery. Although the degradations in the recovered param- 
eters are well above the corresponding quoted statistical 
errors, we do not find a particularly critical effect on the 
various parameters, except for the parameter r— 1, affected 
by an error of 10% by using Jupiter, for the adopted level 
of pointing uncertainty. On the other hand, the analysis 
of the x 2 /DOF (~ 1.3) clearly flags this problem. In this 
circumstance, it may be then particularly safe to take ad- 
vantage of the Saturn transits, that, although intrinsically 
noisier, are significantly less affected by this kind of sys- 
tematics. We can then exploit the Jupiter transits to check 
the fit results where the signal to noise ratio from Saturn 
transits decreases. 

The reported values of x 2 /DOF can be understood 
on the basis of simple considerations. The derivative of a 
Gaussian response, J = exp[— (9 /a) 2 / 2], with respect to 
the angle 9 from the beam center at 9 ~ a is given by 
dJ I J ~ —dO /a; for Jupiter (Saturn), the signal at 9 ~ a 
is - 22 mK (~ 5 mK) for a beam with a FWHM - 33'. 
Then, a pointing error d9 ~ 1' about 9 ~ a implies a 
signal change of ~ 1.6 mK (~ 0.35 mK ). Let's assume 
such value of error in a squared region about the beam 
center with a 2a side, containing ~ 32 samplings in the 
current simulations. The fit uses ~ 10 4 samplings with 
a white noise sensitivity of ~ 0.17 mK. We find then a 
X 2 /DOF of (10 4 x 0.17 2 + 32 x 1.6 2 )/(10 4 x 0.17 2 ) - 1.3 
for Jupiter (and analogously ~ 1.01 for Saturn), in good 
agreement with the values reported in Table 2. 

3.1.2. Spacecraft rotation 

We have also verified the effect the spacecraft rotation dur- 
ing single sample integration. Strictly speaking, the signal 
inside each sampling time changes according to the beam 
response. We can take this effect into account by oversam- 
pling and then averaging the simulated data inside each 
sampling time and by equivalently implementing in the fit 



procedure the recovery of the intrinsic (i.e. non affected 
by rotation) optical beam parameter. We have here imple- 
mented the impact of the satellite rotation only in the data 
flow generation, in order to estimate the error introduced 
in the beam parameter recovery by a simple fit procedure 
that does not take into account this effect. The results are 
again reported in Table 2 for the first transit of Jupiter 
and Saturn in the case of an elliptical beam. As expected, 
the effect on beam position and inclination is negligible. 
Since in the present test the minor axis is taken along 
the scan circle direction, the beam is reconstructed with 
a smaller ellipticity and a larger FWHM with respect to 
the input values, with no indication of a worsening of the 
X 2 /DOF. Clearly, in the general case, the final effect will 
depend on a combination of beam ellipticity and inclina- 
tion with respect to the sky scan direction. The beam pa- 
rameters recovered by neglecting the spacecraft spinning 
within the sampling time in the fit procedure may also 
be seen as characterizing the "effective" properties of the 
beam in presence of spacecraft rotation. Again, the "effec- 
tive" values of a reported in Table 2 are in agreement with 
simple analytical considerations. In fact, the effective reso- 
lution which takes into account the intrinsic beam optical 
resolution and the beam smearing introduced by the satel- 
lite rotation is well approximated by er e / / = \/ a 2 + 9 2 /12 
where 9 S is the angle (~ 5' for the present simulations) in 
the sky described by the beam axis during the sampling 
time. 

3.2. Recovery of the main beam shape 

We consider here the capability to reconstruct in flight the 
detailed shape of the main beam by using Jupiter transits. 
We apply our method to a beam @ 30 GHz simulated 
through the GRASP8 code for the Alcatel easel telescope 
configuration (see Fig. 3), an aplanatic solution like those 
suggested by Mandolesi et al. 2000b in order to minimize 
the coma distortion and render the main beam shapes 
close to ellipses. In this case we adopt a simple scanning 
strategy with a = 80°, the same scan angle considered in 
the Alcatel easel design, shift the spin axis of 2.5' every 
hour and oversample each scan circle every ~ 11' (i.e. 
rougly 3 points per FWHM @ 30GHz). 

By inverting Eq. ([!]) we derive the beam pattern shape. 
Fig. 4 show the results obtained @ 30 GHz by consider- 
ing the sensitivity of the two radiometers coupled to a 
single beam and a single Jupiter transit: the main beam 
shape can be directly recovered down to ~ —20 dB with 
good accuracy. The result clearly improves by adding three 
Jupiter transits and taking into account the possibility 
of averaging the recovered shapes of the two equivalent 
beams at the same frequency, as shown in Fig. 5. 

By fitting the TOD with the method described in 
Sect. 3.1 we can derive the beam parameters under the 
assumption of a circular or bivariate Gaussian shape. The 
results are reported in Table 2. We show also respectively 
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Table 1. Input parameters of symmetric and elliptical beam for the considered planet transits. 



Input values 
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Xt • 10' 


Vt ■ 10" 


Event 


(deg) 


(arcmin) 






(mK) 








circ. 


ellipt. 




circ. 


ellipt. 








Jupiter (I) 












35.8297 






Jupiter (II) 







14.01381 


1 


1.3 


35.4743 


5.76035 


7.91971 


Saturn (I) 






7.8543 


Saturn (II) 












7.8229 







* The value rk depends on the distance between the spacecraft and the planet which slightly varies 
between different pointing events. 
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Fig. 3. Contour plot (in dB) of one of the two beams @ Fig. 4. Reconstruction of the beam of Fig. 3 by using a 
30 GHz as simulated through the GRASP8 code for the single Jupiter transit and the sensitivity of the two ra- 



Alcatel easel telescope configuration. 



diomcters coupled to the beam. 



in Fig.s 6 and 7 the comparison between the input beam 
shape and the recovered one for the case of circular and 
bivariate Gaussian approximation. As evident, the bivari- 
ate approximation represents a significant improvement 
with respect to the circular one, as indicated by the value 
X 2 /DOF, respectively 1.232 and 8.094. The agreement 
between the bivariate approximation and the simulated 
beam results quite good down to ~ —25 dB. At lower re- 
sponse levels the Gaussian behavior begins to significantly 
underestimate the beam response, even in the bivariate 
approximation. 

Finally, we simulate Jupiter transits as seen by an el- 
liptical beam with the same parameters derived from the 
fit of the simulated beam, but in the case of the current 
simple baseline scanning strategy defined by a = 85° , spin 
axis shift of 2.5' every hour and by approximately 3 sam- 
plings per FWHM (i.e. one point every ~ 11' along a scan 
circle @ 30 GHz). As expected, the beam parameters are 



recovered with an accuracy essentially unchanged with re- 
spect to the results of Sect. 3.1 (sec the last row of Ta- 
ble 2), in spite of the less regular grid in the xt — Ut plane, 
less refined (by a factor ~ 4) along the scan direction than 
along that of the spin axis repointing. 

3.3. Circular versus elliptical approximation of the beam 
shape 

A crucial point for CMB anisotropy experiments is to es- 
timate the impact of the quality of the main beam recon- 
struction on the data analysis and, ultimately, on the sky 
maps that can be obtained and on the related science. An 
accurate discussion of this problem is not the aim of this 
work. Nevertheless, we want address here this argument. 
By including also the Galaxy emission, only at low Galac- 
tic latitudes and at the lowest LFI frequencies the impact 
of elliptical main beam distortions significantly increases 
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Table 2. Recovery of the beam parameters from the considered planet transits in presence of pure white noise and 
taking also into account two kinds of systematical effects. Circular, elliptical and simulated beams as well as different 
values of the boresight angle a are considered. 



Event 


e 

(deg) 


a 

(arcmin) 


r 


(mK) 


x* T ■ 10 2 


Vt ■ 10" 


X7DOF 


Circular beam 


Jupiter (I) 




14.0071 
±0.0089 


1.0039 
±0.0012 


35.877 
±0.032 


-5.7602 
±0.0004 


7.9192 
±0.0004 


0.994 


Jupiter (II) 




14.0126 
±0.0092 


1.0010 
±0.0013 


35.486 
±0.033 


-5.7603 
±0.0004 


7.9200 
±0.0004 


0.997 


Saturn (I) 




13.996 
±0.042 


1.0067 
±0.0061 


7.839 
±0.033 


-5.7659 
±0.0020 


7.9204 
±0.0015 


0.998 


Saturn (II) 




14.034 
±0.042 


1.0078 
±0.0061 


7.797 
±0.033 


-5.7637 
±0.0020 


7.9209 
±0.0015 


1.007 


Elliptical beam 


Jupiter (I) 


-0.0059 
±0.0024 


14.0095 
±0.0088 


1.3023 
±0.0016 


35.876 
±0.032 


-5.7602 
±0.0004 


7.9193 
±0.0003 


0.995 


Jupiter (II) 


0.0011 
±0.0025 


14.0102 
±0.0092 


1.2993 
±0.0017 


35.491 
±0.033 


-5.7603 
±0.0004 


7.9199 
±0.0003 


0.997 


Saturn (I) 


-0.007 
±0.011 


14.004 
±0.042 


1.3055 
±0.0078 


7.837 
±0.033 


-5.7661 
±0.0020 


7.9199 
±0.0015 


0.999 


Saturn (II) 


0.009 
±0.011 


14.041 
±0.042 


1.3073 
±0.0078 


7.796 
±0.033 


-5.7645 
±0.0020 


7.9207 
±0.0015 


1.007 


Elliptical beam: effect of 1' pointing error 


Jupiter (I) 


0.0288 
±0.0025 


13.7895 
±0.0089 


1.3424 
±0.0017 


35.884 
±0.032 


-5.7576 
±0.0004 


7.9227 
±0.0003 


1.329 


Saturn (I) 


-0.005 
±0.011 


14.004 
±0.043 


1.3060 
±0.0078 


7.819 
±0.033 


-5.7589 
±0.0020 


7.9118 
±0.0015 


1.011 


Elliptical beam: effect of neglecting spacecraft motion 


Jupiter (I) 


0.00367 
±0.0025 


14.0697 
±0.0089 


1.2904 
±0.0016 


35.584 
±0.032 


-5.7602 
±0.0004 


7.9198 
±0.0003 


0.999 


Saturn (I) 


-0.007 
±0.012 


14.053 
±0.043 


1.2964 
±0.0078 


7.782 
±0.033 


-5.7661 
±0.0020 


7.9200 
±0.0015 


0.999 


Alcatel easel beam: fit with circular beam - realistic simple scanning strategy, a = 80° 


Jupiter (I) 




14.9583 
±0.0089 




31.106 
±0.026 


-5.5751 
±0.0003 


7.2991 
±0.0004 


8.094 


Alcatel easel beam: fit with elliptical beam - realistic simple scanning strategy, a — 80° 


Jupiter (I) 


20.97 
±0.11 


15.0102 
±0.0087 


1.3602 
±0.0016 


31.360 
±0.026 


-5.5756 
±0.0003 


7.3001 
±0.0004 


1.232 


Above elliptical beam: fit with elliptical beam - baseline simple scanning strategy, a = 85° 


Jupiter (I) 


21.01 
±0.10 


15.0209 
±0.0083 


1.3590 
±0.0015 


32.241 
±0.025 


-5.5792 
±0.0003 


7.2794 
±0.0004 


0.997 



(by a factor ~ 3) with respect to the case of a pure CMB 
fluctuation sky (Burigana et al. 1998) whereas only very 
small effects are added by the combined effect of realistic 
main beam distortions and extragalactic source fluctua- 
tions (Burigana et al. 2000a). In addition, the main beam 
distortion impact moderately increases with the FWHM 
(Burigana et al. 1998) for the LFI resolution range. There- 
fore, we will consider here for simplicity a pure CMB fluc- 
tuation sky in a reference case of a standard CDM model 
approximately COBE/DMR normalized and consider for 
the present analysis the 30 GHz channel. In Fig. 8 we 
report the difference between the signals obtained by con- 
volving the sky with the bivariate main beam and with 
the circular main beam by assuming in each case the cor- 



responding set of parameters that fits the simulated beam 
(see Table 2). The rms of the temperature differences is 
~ 4.1/iK. A similar comparison but between the circular 
and the simulated beam (see Fig. 9) gives a very similar 
rms value, ~ 4.7^iK. This suggests that the ellipticity is 
the most relevant feature of the main beam shape. In fact, 
the comparison between the convolutions with the bivari- 
ate Gaussian beam and the simulated one (see Fig. 10) 
gives a rms temperature difference of ~ 1.7/xK, a fac- 
tor ~ 3 smaller than that obtained by approximating the 
beam with a symmetric Gaussian shape. 

Of course, when the Planck optical design will be 
settled and the main beam patterns computed through 
optical simulation codes it would be possible to search for 
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Fig. 5. Reconstruction of the beam of Fig. 3 by using three Fig. 7. Comparison between the simulated beam and its 

Jupiter transits and two radiometers per beam and by approximation in terms of elliptical Gaussian beam (see 

taking advantage from the identical optical properties of also Table 2). 
the two 30 GHz beams. 
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Fig. 8. Difference between the TOD from the elliptical 
Fig. 6. Comparison between the simulated beam and its beam and thc drcular beam that fit thc input beam for 

approximation in terms of circular Gaussian beam (see tij ree scan circles 
also Table 2). 



analytical descriptions of main beam shapes that might 
improve the bivariate Gaussian approximation. Anyway, 
we have proved here that this representation allow to reach 
the ~ /xK accuracy level in the treatment of Planck-LFI 
TOD's. 



4. Straylight from Moon, Earth and Sun 

By using the same method described in Sect. 2 and the full 
antenna pattern shown in Fig. 1 (de Maagt et al. 1998), we 
are able to quantify the impact of the internal Solar Sys- 
tem objects to the Planck-LFI observations. The Moon, 
the Earth and the Sun are the only internal bodies that 
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Fig. 9. Difference between the TOD from the circular 
beam and the simulated beam for the same three scan 
circles of Fig. 8. 
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Fig. 10. Difference between the TOD from the elliptical 
beam and the simulated beam for the same three scan 
circles of Fig. 8. By comparing this figure with Fig. 9 the 
improvement with respect to the circular approximation 
of the main beam is evident. 

may introduce appreciable straylight contamination. The 
solid angle of these objects, although not negligible, is very 
small compared to the angular scale for which significant 
variations of the far pattern response occur; this make the 
use of Eq. (2) accurate enough for the present purposes. 



Of course, the level of straylight contamination from 
these bodies depends also on the effective scanning strat- 
egy. We studied a Lissajous orbit around L2 for the 
Planck surveyor (Bersanelli et al. 1996) by considering 
case i) the simple scanning strategy with the spin axis al- 
ways parallel to the Sun-spacecraft direction and case ii) 
a precession motion of the satellite spin axis like that dis- 
cussed in Sect. 2.1. [For this straylight analysis we assume 
a = 80° , the value adopted in the optical computation by 
de Maagt et al. 1998]. 

By coadding the TOD computed as described in Sect. 2 
we have produced nearly full sky maps of averaged stray- 
light contamination (sec Burigana et al. 2000b for fur- 
ther details). [We adopt here the HEALPix pixelisation 
by Gorski et al. 1998]. We find that in any case this ef- 
fect is very small, with the maximum contamination level 
always below 0.15 /iK. Although small, this effect has to 
be considered with caution. In fact, this straylight con- 
tamination is produced by the very low response antenna 
in the "antispillover" region, where an accurate antenna 
response measurement is extremely difficult. 

5. Discussion and conclusions 

We implemented the Planck flight simulator (Burigana 
et al. 1997, 1998, Maino et al. 1999) to properly discuss the 
impact of the Solar System main bodies on the Planck 
observations. In particular, we focused on the problem 
of the in-flight reconstruction of the main beam of the 
Planck-LFI antenna patterns. To do so, we simulate in 
details the transits of Jupiter and Saturn in the field of 
view of the Planck-LFI, 30 GHz beam. The method can 
be easily extended to the other Planck channels. Our 
analysis shows that, using Jupiter, we can recover in flight 
the main beam response down to ~ —(25-^32.5) dB, where 
the signal to noise ratio approaches unity. 

Both symmetric and non symmetric beams have been 
considered; in the latter case we assumed a simple el- 
liptical shape, as suggested by recent optical simulations 
(Mandolesi et al. 2000b), but the method can be general- 
ized to more refined parametrizations. 

We have demonstrated that the key parameters of the 
main beam (resolution, cllipticity, position and inclination 
on the plane of Planck field of view) can be simultane- 
ously recovered with high precision by fitting the planet 
transit signal. Of course, the larger signal to noise ratio 
of Jupiter (compared to that of Saturn) translates in a 
better parameter recovery, by a factor 3 -j- 5. 

After having considered the idealized case of pure 
white noise, we discussed possible degradations to the 
main beam reconstruction introduced by the effects of 
spacecraft pointing relevant in this context. Spacecraft ro- 
tation does not affect significantly the quality of the re- 
construction; on the contrary, bu neglecting the spacecraft 
rotation in the fit procedure we derive the effective beam 
parameters including rotation smearing effect. The most 
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relevant source of contamination is represented by the spin 
axis pointing uncertainty which can degrade the recovery 
of the beam ellipticity parameter r — 1 by ~ 10% in the 
case of V pointing error when Jupiter is used. In these 
circumstance, clearly recognizable through the increase of 
the x 2 / DOF, we find very advantageous to consider Sat- 
urn, that produce information less sensitive to this kind 
of systematics. 

To complete the analysis, we considered also the full 
recovery of a simulated beam shape for an aplanatic con- 
figuration of the optical design recently studied by Alca- 
tel. The beam shape can be reconstructed with good accu- 
racy and resolution down to —25 dB. Of course, somewhat 
lower response levels at larger angles from the beam centre 
can be measured in flight by relaxing the requirement of 
recovering the beam shape with the same resolution nec- 
essary to accurately describe the more central regions of 
the main beam. 

We addressed the study of the impact of the quality 
of the main beam reconstruction on the treatment of LFI 
TOD's. We shown that a bivariate Gaussian approxima- 
tion represents a significant improvement with respect to 
a symmetric representation. A rms difference of ~ 1.7/^K 
is found between the TOD obtained by convolving a CMB 
fluctuation sky with a simulated beam or its bivariate 
Gaussian approximation as derived by the in-flight recon- 
struction, a result better by a factor three with respect to 
the case of the symmetric approximation. 

The possibility to combine a very accurate in-flight cal- 
ibration by using the CMB dipole (Bersanelli et al. 1997) 
and the good accuracy in the recovery of the maximum 
signal (the parameter in Table 2) at the planet transit, 
offers a good chance of measuring the intrinsic planet tem- 
peratures at millimctric wavelengths with an accuracy at 
% level, the main source of error being the uncertainty on 
the integrated antenna pattern response. This represents 
an interesting byproduct of Planck observations. 

Finally, we have shown that the effect of Sun, Earth 
and Moon in the far sidelobes produces sub-/xK effects 
independent of the details of the scanning strategy. 

To summarize, at least at 30 GHz, observation of exter- 
nal planets offers an accurate, robust and simple method 
to reconstruct in flight the main beam properties under 
very general conditions. 



LFI27 for the Alcatel easel configuration with the high 
resolution grid adopted in this study. 
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